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Abstract 



We consider the effects on phonon dynamics of spin-lattice coupling within 
the slave-boson mean-field treatment of the extended t-J model. With no ad- 
ditional assumptions the theory is found to give a semi-quantitative account 
of the frequency and linewidth anomalies observed by Raman and neutron 
scattering for the 340cm^^ Big phonon mode in Y Ba2Cu^0'j at the super- 
conducting transition. We discuss the applicability of the model to phonon 
modes of different symmetries, and report a connection to spin-gap features 
observed in underdoped YBCO. The results suggest the possibility of a unified 

understanding of the anomalies in transport, magnetic and lattice properties. 
PACS numbers: 74.20.Mn, 74.25.Ha, 74.25.Kc 
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Anomalies in the transport and magnetic properties of high-temperature superconduc- 
tors, observed in both normal and superconducting states, have attracted much interest [|I| , 
since these phenomena are considered to be manifestations of a metallic state arising near the 
Mott transition due to strong correlations. Although there is as yet no general consensus on 
the appropriate theoretical description of this anomalous metallic state, the approach based 
on the t-J model |^ treated by the slave-boson mean-field theory has given important 
clues to the understanding of both transport and spin excitation properties . These 
include the temperature-dependences of the resistivity and Hall coefficient, and the different 
temperature-dependences P,p^ of the shift and rate of nuclear magnetic resonance between 
high- and low-doping regions, especially the "spin-gap" phenomenon 1^,^] first noted by Ya- 
suoka 0. There have been various experimental reports |[TlHl3| which indicate anomalous 



temperature-dependences also in the frequency shift and linewidth of phonons around the 



superconducting transition. Further, studies by neutron scattering [|14|,|T5|, EXAFS [[16 
ion-channeling and ultrasound measurements ||18| all provide strong evidence of a link 
between lattice anomalies and superconductivity. Although some aspects of the phonon 
problem have been studied theoretically |T^, to our knowledge there have not so far been 
any efforts to understand all of these interesting low-lying excitations on a unified basis. 

We consider the spin-phonon coupling arising naturally within the extended t-J model of 
a single Cu02 layer, which has been shown at the mean-field level to give a good account 
of many features of the spin excitations in both La2-xSrxCu04 (LSCO) and YBa2Cu30'i_s 
(YBCO). The Hamiltonian is 

H = tijCilgCijs + JijSi.Sj, (1) 

ij (ij) 

where the Hilbert space is that without double occupancy, tij corresponds to the transfer 
integrals used to reproduce the Fermi surface, and the superexchange interaction Jij, which 
is assumed to be finite only between nearest neighbours, has been taken as a constant in 
previous treatments. 

In YBCO, the Cu02 layer is "buckled", by which is meant that the oxygen atoms 0(2) 



and 0(3) lie out of the plane of the Cu atoms, as shown schematically in Fig. 1. Then tij 
and Jij between the nearest-neighbor Cu sites have contributions linear in the magnitude of 
the oxygen displacement along the c-axis, uf {i refers to the Cu sites on the square lattice, 
and a to either 0(2) (a = x) or 0(3) (a = y)), given by Uj = t[l — \t{uf/a)], where 
a and At are the distance between Cu sites and the coupling constant, respectively, and 
similarly for Jij with coupling Aj. A microscopic estimate of At and Aj requires several 
steps. From the structural data of Ref. 12^], the equilibrium O displacement is uq = 0.256A 



(neglecting henceforth the 5% anisotropy between x and y) ; although the degree of buckling 
is small, uo/a <^ 1, its inclusion is crucial in providing a coupling which is strong and linear. 



Following Ref. ||2T|, the dependence on interatomic separation of the transfer integral to- {tn) 



between Cu-.d^^-y^ and 0:pa (O-.p-^) is given by 
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where ^Vpda {Vpdn) is the transfer integral between dx2_y2 (d^z) and p^- (Pz) orbitals with 
separation d = [(|a)^ + 'U^] along the x-direction. Taking Vpd{d) oc c?^^/^ and writing 



u = Uq + 6u, where 6u is the oscillation amplitude, gives t„{u)/t„{uo) = 1 — 2.035^/0 
and tT,{u) /taiuo) = 1.536u/a. The final requirement is the dependence on and t-,, of tij 
and Jij. For Jij, the perturbative expression J = 4 (t^ — 2t^t^) / A U (cf. |2^) would give 
a large coupling constant A = 10.4. However, use of the lowest-order perturbation form 
is not well justified, and investigations of the influence of higher-order terms find for the 
parameters of the Cu02 layer the effective relationships J oc tp^ with x ^ 2.3 [^, and 
t oc with ?/ ~ 1.0 |2^; the result for y requires consideration of the many inter-site 
transfer integrals which contribute to the hopping of a Zhang- Rice singlet. The exact values 
of X and y, and particularly of the powers y' and y" corresponding to the extended hopping 
terms t' and t", are quite sensitive to both the parameter choice and to the symmetry of 
the phonon mode EM. This treatment disregards changes in A, which we believe to be 



appropriate for the local, screened deformation processes under investigation. By contrast. 



a (i-dependence of A is required to account for the weak relation J cc d ",4<a;<6, found 
by pressure- and substitution- induced variation of the bond length where in addition 



the Madelung energy is altered |2^. Mindful of these uncertainties, we proceed by taking 
the values x = 2.0, y = 1.0 so that Aj = 5.2 and At = 2.6 = |Aj, and neglect modulation 
of the extended tij terms for the purposes of the current study. These estimates have also 
been found to give good agreement with the measured isotope shift in YBCO |2^, a result 
(deferred to a future publication) which provides an independent indication of their validity. 

The terms in H (Eq. (|l])) describing the coupling of uf to the spin degrees of freedom 
may be rewritten in the slave-boson, mean-field approximation as 

-EE {t)^t{^){b^bla)X^,^+a (4) 
i a=x,y 

+ l>^^^(¥) [{xl+a)Xi,+a + 2{At)A,, + h.C.]}, 

where Xi,j = J2s fLfjs and Aij = [fi^fji — fiifji) / V^- We do not consider the phonon-holon 
coupling vertex which is also contained in the t term, because in the normal state this will 
vanish at g = 0, while in the Bose-condensed state the ho Ion has no dynamics. 

Here we will be concerned mainly with the T-dependence of the frequency shift and 
linewidth of the 340cm^^ Big mode in Y Ba2Cu^0-j, which is an out-of-phase oscillatory 
motion of only the planar oxygen atoms (wf = — "uf), and is also out of phase between the 
planes of the bilayer This mode has attracted experimental interest because it shows 



the largest effects at the superconducting transition, and there are available detailed Raman 
TTHT3[ (g = 0) and inelastic neutron scattering (also g 7^ 0) data. We will also show 



results for Aig- or A2u-symmetric (in-phase, = u\) oscillations, and discuss the relevance 
of the model for these. 

The effect of the coupling on the dynamical properties of the phonon is calculated from 
the lowest-order spinon correction to the phonon self-energy Ilph (Fig. 2). Second-order 
perturbation theory in terms of the coupling results in a frequency shift 5uj =ReI\.ph for 
l^cul ^ ojq, which is given at g = by 
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o 4 ^ 1 tanh fl^) 
6u = c{Xjjfi-yBk o (5) 

where c = (^)%^) = 1.18 x 10~\ {u^) = = (0.055A)2, in which uo = SAOcm-^ is 

taken for the phonon frequency of interest and M is the mass of the O atom, Ek is defined 
by Ek = l^l + A^]^^^, with the spinon band energy relative to the chemical potential and 
Afc = — ^j^JA(cosfc^ — cos ky) the singlet order parameter, and is the form factor 

Bk = 2A^ (7fc^fc + ^xvlf Big, B2u modes 

Bk = 2A^r]l (a + xX7fc)' ^ig, Mu modes, (6) 

with 7fc = coskx + cosfcj^, 77^ = cosfc^; — cosky, A = (Aj^) and x = {Xij) + jj- Oiie observes 
in Eq. @ that 5^ oc A^, and so there is no phonon energy correction due to spin coupling 
in the normal state; this result is a special feature of the u ^ phonon mode we consider at 
g = 0, and of the current level of approximation. The superconductivity- induced correction 
to the linewidth F is computed similarly from ImUph. 

A detailed numerical analysis of the self-consistent mean-field equations has been carried 
out to compute the temperature-dependences of A, x and chemical potential for several 
choices of doping 6, for the transfer integrals t = 4 J, t' = — |t and t" = |t, appropriate 
to YBCO. With these parameters, the frequency-dependence of 6u!, given by Eq. (||), has 
been evaluated numerically for several temperatures, and the results are shown in Fig. 3 
for g = and 6 = 0.2 at T/Tryb = 0.2, where T^yb = 0.069J is the onset temperature 
for the singlet RVB order parameter A. All calculations were performed with an assumed 
Lorentzian broadening of the spinon spectrum F = 0.12kBTRVB, and are found to be quite 
insensitive to the value of this parameter. From the tu-dependent denominator in Eq. (^, 
the frequency where 6uj of the Big mode changes sign is an approximate measure of the 
value of 2Afc(T) near the (vr, 0) points, where the gap is maximal: at higher temperatures 
(but below Trvb) the tu-dependence of 5uj is qualitatively the same, but the sign-change 
occurs at lower frequencies. The experimental mode frequency ujq for the Big mode is near, 
but just below, the low-T crossover, so this mode can be expected to show a maximal effect. 



The linewidth broadening ST has the form to be expected for the corresponding imaginary 
part, namely a peak at the crossover frequency. 

The temperature-dependence of 6uj is shown in Fig. 4(a) for g = and uj/J = 0.25, 
a value in reasonable quantitative agreement with the 340cm~^ mode: while modes with 
frequencies considerably less than 2A have a sharp transition, those close to it show clearly 
that as uj approaches the maximal value of 2Afc(T = 0), the shift in frequency occurs at 
a temperature somewhat below Trvb- This corresponds well to the observations of Ref. 
where the full frequency shift develops over a range of temperatures below the onset. 



Comparison with the experimental result |[T2| for the broadening of the Big mode again 
shows a good agreement in sign and magnitude {5uj ~ 5r ~ O.Olcuo) for a frequency close 
to ujq (Fig. 4(b)), while at lower frequencies ("off resonance") 5V is suppressed. We have 
extended our analysis to the case where the phonon has a finite wavevector q, finding smaller 
but sharper anomalies in further qualitative agreement with experiment; these calculations 
will be presented elsewhere. Quantitatively, the magnitude of the effects given by the model 
with the chosen values of Aj and Aj is within a factor of 1.5 of the Raman measurements Wl 



on Y Ba2Cu^0j] this degree of correspondence, as well as that in ujq above, can be regarded 
as satisfactory within a mean-field treatment using no adjustable parameters. 

In Figs. 3 and 4 are shown not just the Big mode, but also the results for a mode of in- 
phase 0(2) and 0(3) oscillations, which would correspond to Aig (Raman active, 440cm~^) 
or (infrared- active, 307cm~^) symmetry. At the level of the current approximation these 
have negligible anomalies, a qualitative difference from Big which arises from the form factors 
Bk (§) with (i-wave singlet pairing; we note that extended s-pairing gives an immeasurably 
small frequency shift for both types of mode symmetry. However, experimental studies 



of the mode |]TT|, whose frequency is close to that investigated, reveal a strong shift 
5uj/ujq ~ 1%. We believe that the model cannot reproduce this result because it contains 
nothing to account for the charge transfer between the Cu02 planes of a bilayer which 
accompany this mode in YBCO, and so is inapplicable. Developing the theory to include 
interplane charge motion is beyond the basic t-J framework, and will be pursued in a 
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subsequent publication. 

The 193cm~^ mode is an out-of-phase oscillation of planar O (in-phase between the 
planes of the bilayer) involving little charge motion, and thus a single-layer model is expected 
to be valid. Because this mode is i.r.-silent, it may be studied only by neutron scattering, and 
investigations on Y Ba2Cu^0j have recently been performed The observed anomalies 
are sharp and occur close to Tc, with a relative magnitude 5uj/ujq ~ 1% similar to those of 
the Big mode, features which are indeed well described by the current model at this lower 
frequency. Finally, the effects studied should be strongly suppressed in the i?-symmetric 
phonon modes, where O displacements are parallel to the plane and so have little influence 
on the interaction Jjj, a feature also largely in agreement with experiment. 

These results raise the interesting possibility of probing the spin-gap behavior observed 
in the temperature-dependence of the NMR relaxation rate by considering the frequency 
shifts of particular phonons in members of the YBCO class in the low-doping regime, where 
the rate has been found to be maximal at some value To above the superconducting critical 
temperature T^. There exist already several experimental reports [pT|-[TB|] of anomalies in 
the frequency shift well above T^, and near the temperature where the NMR rate exhibits a 
maximum, while recent, highly accurate Raman studies of phonon anomalies in underdoped 



YBCO compounds |^ show clearly the onset of a frequency shift at some Tq ~ I'bQK, 
followed by growth of this shift as temperature is lowered, until a saturation below T^. Such 
features may be understood on the basis of the mean-field phase diagram of the extended 
t-J model 0, in which T^yb is indeed higher than only in the low-doping region, and 
they would be interpreted as the onset of singlet-RVB order around Tq, i.e. its identification 
with the crossover temperature T^vb [Q- This result is consistent with the observation that 
the specific heat anomalies at Tc in these cases are qualitatively different from those near 
optimal doping p2| , p3| . We suggest that the existence of two temperature scales is the most 
likely explanation for the contrasting low-g behavior of the energy shifts in nominal O6.92 
and O7 compounds observed in Ref. 

In summary, we have proposed a theory of spin-phonon coupling which accounts very 



well for the anomalies observed in phonon modes in typical high-Tc cuprates of the YBCO 
class, based on the mean-field approximation to the extended t-J model. We beheve that 
this model is a useful step in constructing a coherent theory of the spin, transport and lattice 
properties in the anomalous metallic state, which has a strong bearing on high-temperature 
superconductivity. 
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FIGURE CAPTIONS 

Fig. 1: Schematic representation of Cu02 layer, showing "buckhng" deformation of the 
equihbrium positions of 0(2) and 0(3) atoms out of the plane of the Cu atoms, appropriate 
for YBCO. a, b and c represent crystal axes. 

Fig. 2: Diagrammatic representation of the lowest-order contribution to the phonon self- 
energy Uph due to coupling to spinous. 

Fig. 3: Phonon frequency shift 5u; for (o) and (x) modes at 5 = 0, as a func- 
tion of frequency at T = 0.2TRyB. The arrow indicates the frequency whose temperature- 
dependence is illustrated in Fig. 4. 

Fig. 4: Phonon frequency shift Slu (a) and linewidth broadening ST (b) for Big (o) and 
(x) modes at g = 0, as functions of temperature for mode frequency cuq — 0.25 J. 
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